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We investigate the fractionalized Skyrmion excitations induced by spin-orbit coupling in rotating 
and rapidly quenched spin-1 Bose-Einstein condensates. Our results show that the fractionalized 
Skyrmion excitation depends on the combination of spin-orbit coupling and rotation, and it orig- 
inates from a dipole structure of spin which is always embedded in three vortices constructed by 
each condensate component respectively. When spin-orbit coupling is larger than a critical value, 
the fractionalized Skyrmions encircle the center with one or several circles to form a radial lattice, 
which occurs even in the strong ferromagnetic/antiferromagnetic condensates. We can use both the 
spin-orbit coupling and the rotation to adjust the radial lattice. The realization and the detection 
of the fractionalized Skyrmions are compatible with current experimental technology. 
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Introduction. — Spin-orbit coupling (SOC) describes 
the interaction of a particle's spin with its motion. The 
particular form of SOC can be of either Rashba [l[ or 
Dresselhaus type. SOC in an electronic system [sl 
is able to serve as a spin filter or a Stern- Gerlach ap- 
paratus. And it is crucial^or the spin-Hall effect 0, IH 



and topological insulators [6|-l9| . Recently, spin-orbit cou- 
pled Bose-Einstein condensate (BEC) has been realized 
in NIST's experiment [lo| for the first time. Unlike pre- 
vious experiment, their work has factually explored the 
bosons system in the non-Abelian gauge field EMI. 
This opens up a new avenue in cold atom physics and 
attracts many attentions. 

Motivated by NIST's experiment [lo|, several recent in- 
vestigations about bosons with SOC have presented some 
nontrivial new structures such as stripe phase [l8l-[2l| and 
half-quantum vortex state 2qI-[22| . Specially, the combi- 



nation effect of SOC and rotation on pseudo spin-| BEC 
has been shown to be able to generate various vortex 
structures [H, [1^ . These impressive results enrich the 
phase diagram of BEC system. However, how can SOC as 
a new effect on spinor BEC produce previously unknown 
types of topological excitations such as new Skyrmion? 

In this Letter, we explore how SOC induces the frac- 
tionalized Skyrmions in rotating spin-1 BEC. We find 
that the fractionalized Skyrmion excitation which is re- 
lated to a three-vortex structure exists in the rotating 
spin-1 BEC with SOC. Unlike previous Skyrmion-like ex- 
citations, the generation of the fractionalized Skyrmion 
excitation must depend on the combination of SOC 
and rotation. In the absence of rotation no fractional- 
ized Skyrmion emerges, and this system can degener- 
ate into the stripe phase. Without SOC, the integral 
Skyrmion excitations can appear in the ferromagnetic 
BECs. When SOC is larger than a critical value, the 
fractionalized Skyrmions encircle the center of the sys- 



tem with one or several circles to form a radial lattice. 
These phenomena can occur even in the strong ferromag- 
netic/antiferromagnetic BECs. We can adjust the frac- 
tionalized Skyrmions lattice by changing the strength of 
SOC as well as rotation. 

Model and equation. — Considered a quenching process, 
the dynamics of a spin-1 BEC can be described by the 
stochastic projected Gross-Pitaevskii equation |25|422|- 
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where \E'j(j = 0, ±1) denotes the macroscopic wave func- 
tion of the atoms condensated in the spin state \F = 
l,mF = j), and H,^, = [-^ + V{r) + 5„|*n^ . + 
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spin-1 matrice, Q is the rotation frequency, Lz [L^ = 
—ih{xdy — ydx)] is the z component of the orbital an- 
gular momentum, and tz denotes the strength of SOC. 
Meanwhile, T is the final temperature, /i is the chemical 
potential, is the growth rate for the jth component, 
and dWj is the complex Gaussian noise. The projection 
operator V is used to restrict the dynamics of the spinor 
BEC in the coherent region. 

In numerical simulations, the initial state of each 
is generated by sampling the grand canonical ensemble 
for a free ideal Bose gas with the temperature Tq and 
the chemical potential /ij^o- Meanwhile, the condensate 

band must lie below the energy cutoff Er > Ek = — 2m * 
Noting, k = 27r{nx^ny)/L, where rix, riy are integers and 
L is the size of the computation domain. To simulate the 
quenching process, the final temperature and the chem- 
ical potential of the noncondensate band are altered to 
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FIG. 1: (color online). The densities and phases for the spinor 
BEG of ^^Rb with SOG k when the system reaches the equi- 
librium state, (a) k, = 0.1; (b) k, = 0.2; (c) k = 0.5; (d) 
K — 0.7; (e) HI — 1.0. Here, Q = O.Scj, ao = 101. Sas and 
a2 = 100. 4aB. Noting, the fifth and sixth columns are the 
phases of = — 1 and = 1 components, respectively. 



the new values T < Tq and /i > /ij^o- Furthermore, we 
use the oscillator unit in the numerical computations. 

The length, time and energy are scaled in units of a/ — , 

cj~^ and hw^ respectively. In the simulations, the total 
number of the modes are fix •, 

= 300, the energy cutoff 
is chosen at Uxc^ ^lyc = 150, the final temperature T is 
lOni^^, the chemical potential ji = 25fia;, and = 0.03. 

Fractionalized Skyrmion lattice in rotating spin-1 BEC 
with spin- orbit coupling. — We begin with the spinor BEC 
of ^^Rb [28], which is ferromagnetic (FM) {qs < 0). Fig- 
ure 1 displays the densities and phases obtained under 
various strength of SOC. For a very weak SOC = 0.1), 
the patterns are irregular. When is over 0.1, the pat- 
terns are relatively regular. The center of the systems 
displays various geometrical structures such as triangle, 
square, heptagon etc. These pictures factually display 
the transition of the patterns as the strength of SOC in- 
creases. Thus, SOC can be used to adjust the pattern in 
the rotating spin-1 BEC. 

Just as previous experiments about the rotating BECs 



251 . |26| , there are some vortices in the three components 
respectively. The fifth and sixth columns indicate the 
phases of = —1 and mi? = 1 components respec- 
tively. Like the vortex lattice in single-component BEC, 
there are some lines where the phases change discontinu- 
ously from red to blue, which corresponds to the branch 
cuts between the phases — tt and tt. The ends represent 
phase defects. All the lines extend to the outskirts of 
the BEC where the density of the BEC is almost negligi- 



ble, and end with another defect which offers neither the 
energy nor the angular momentum to the system. Fur- 
thermore, we also find some peaks accompanying the vor- 
tices, regularly arraying to be triangle, square, heptagon 
e^c, especially in the center of = — 1 component. 

Unlike the periodic vortex lattice [25] or the vortices 
trimers [26], the vortices encircle the center with several 
circles. The number of vortices is 1 or in the cen- 
ter, and it increases as the radius increases. Certainly, 
this phenomenon is not obvious when SOC is very weak 
{hi = 0.1). The fourth column shows the total density of 
BECs. Here, we can distinguish some local minimum of 
densities, especially when tz approaches 1. 

To explore the possible Skyrmion structure, we further 
examine the spin texture 29|, l30| which is parallel to the 
local magnetic moment. The spin texture is defined by 
= E^,n=o,±i^m(^c.)m,n^n/|^P = x,y,z). Fig- 
ure 2(a) shows the spin texture in Fig. 1(c). Here, the 
arrows form a Skyrmion-like circle only in the center of 
the system, namely the black pane region. By calculating 
the topological charge Q = ^ / S • (§| x ^)dxdy, we 
find Q = 1/2 for the structure. Meanwhile, the arrows 
form big rings, whose main direction is marked with the 
blue arrows. But some mutations occur obviously, that 
is, the direction of arrow changes suddenly and the ar- 
rows form a small half circle locally. Additionally, the 
orientations of the arrows suffer a 180° reversal along the 
radius. For clarity, we use a blue pane to point out the 
special structure, where the topological charge is 0.25. 
There are many of these structures rounding the center 
with several circles, i.e., radially arranging in the system. 
We know the spin vector sweeps a whole unit sphere for 
a Skyrmion excitation with topological charge Q = 1; a 
half Skyrmion has the topological charge Q = 1/2, where 
the spin vector covers only half of a unit sphere (3l| . 
Here, we find that the S vector in the black pane covers 
half of a unit sphere, and the S vector in the blue pane 
covers quarter of a unit sphere. Thus, we can call the two 
structures as half Skyrmion and fractionalized Skyrmion 
respectively. How do the half Skyrmion and the fraction- 
alized Skyrmions emerge in this system? Why does it 
only cause one half Skyrmion? 

Figure 2(b) illuminates the relationship between the 
vortices and the fractionalized Skyrmions. The green, 
blue and red spots are the center of vortices formed by 
mF = — 1, ^F = and mi? = +1 components, respec- 
tively. Except for the central area of the BECs, we find 
the positions of vortices in the three components are far 
away from the center with certain order: green, blue and 
red. Meanwhile, the sequence distributes in the whole 
system. Thus, we can view the three-vortex structure 
as a cell. Undoubtedly, the number of vortices in the 
three components approaches 1:1:1. Furthermore, 
each fractionalized Skyrmion accompanies a three- vortex 
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structure. Usually, the Skyrmion-like excitations 
are related to the underlying vortex configuration such 
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FIG. 2: (color online), (a) Spin texture of spinor BEC of 
^^Rb with hv = 0.5 and Q = O.Scj. The color of each arrow 
indicates the magnitude of Sz- The black pane points out a 
half Skyrmion, and the blue pane indicates a fractionalized 
Skyrmion. The blue arrows show the main direction of the 
spin texture, (b) The position of vortices and the spin tex- 
ture. The green, blue and red spots are the center of vortices 
formed by the ttif = —1, tuf = and ttif = +1 compo- 
nents, respectively, (c) A scheme of three-vortex structure. 
The green, blue and red surfaces denote the densities of the 
rriF = TTiF — and ttif = +1 components, respectively. 
The cyan is the total density of the BECs. (d) The position 
of vortices and the spin texture of spinor BEC of ^^Rb with 
K = 0.1 and Q = 0.5uj. Noting, we only mark the vortices 
in y > region in order to illuminate the spin texture and 
position of vortices clearly, (e) The total topological charge 
of system under the equilibrium state as a function of SOC 
strength where Q = 0.5uj. 



as Mermin-Ho [32| and Anderson- Toulouse [33| vortices. 
Here, the fractionalized Skyrmion is related to the three- 
vortex structure. 

Figure 2(c) further indicates a scheme of the three- 
vortex structure with the density distribution. Here, the 
green, blue and red surfaces represent the densities of the 
mi? = — 1, ^F = and = +1 components, respec- 
tively. The cyan represents the total density. There is a 
local density minimum at the position of vortex formed 
by rriF = component when we examine the total den- 
sity. This point is different from the normal coreless vor- 
tex [29I, lo, 32, 33I where the total density has no singu- 



FIG. 3: (color online). The effect of rotation frequency Q 
for spinor BEC of ^^Na with k = 1, ao = 50a b and a2 — 
55as. (a) Q = 0; (b) Q = 0.2cj; (c) Q = 0.5cj. The fourth 
column shows the corresponding spin textures and position of 
vortices. The meanings of the spots and the colored arrows 
are the same as those in Fig. 2(b). Noting, we only mark the 
vortices in ^ > region. 



larity. The = +1 component forms a obvious hump 
at the vortex region of = —1 component, and vice 
versa. Those properties cause a dipole of spin. Generally, 
the dipole is embedded in the three-vortex structure. 

In the absence of SOC, the periodic Skyrmion lat- 
tice can be created in the rotating spin-1 BEC of ^^Rb 
[2^. Figure 2(d) shows the spin textures and the posi- 
tion of vortices when the strength of SOC is 0.1. Here, 
the vortices hardly form the three- vortex structure, es- 
pecially near the center. Additionally, the fractionalized 
Skyrmion lattice is not very obvious. These results fur- 
ther prove that the fractionalized Skyrmion is related to 
the three-vortex structure. The distribution of fraction- 
alized Skyrmion depends on that of vortices. Thus, it is 
not inconceivable that the different strength of SOC will 
cause various fractionalized Skyrmion lattices. 

Figure 2(e) shows that the total topological charge of 
the system varies as SOC is increasing. Because SOC 
changes the integral Skyrmions into the fractionalized 
ones, the total topological charge of FM BECs tends to 
decrease firstly. Then, the value generally gets greater. 
Furthermore, we find the fractionalized Skyrmion lattice 
can also occur in the antiferromagnetic (AFM) BEC, 
where ^5 > [see Fig. 3(c)]. There is no Skyrmion 
excitation appearing in rotating AFM BEC of ^^Na [2^ 
when K ■= (the total topological charge is in Fig. 
2(e)). Added SOC, the total topological charge in- 
creases. This means that SOC enhances the fractional- 
ized Skyrmion excitations. Especially, to obtain the frac- 
tionalized Skyrmion lattice, the strength of SOC must 
exceed a critical value. Here, this value approaches 0.2. 



The effect of the rotation frequency. — Here, we use the 
spinor BEC of ^^Na to illustrate the effect of the rota- 
tion frequency. We only change the rotation frequency 
Q and fix all other parameters to perform the numeri- 
cal experiments. Figure 3 shows the density distribution 
and spin texture under various rotation frequencies. In 
the absence of rotation = 0), there is no vortex ap- 
pearing at all. Each component of the BECs is split 
to be several parallel parts. In fact, these properties 
agree with the stripe phase [l8|-[Il|- The spin texture 
indicates no fractionalized Skyrmion excitations in this 
system. The color straps factually are the low density 
domain of BECs. Added a weak rotation (0.2cj), the 
splitting parts bend and break, and several vortices and 
the three-vortex structure occur. When the rotation be- 
comes faster (0.5a;), the vortex lattice emerges and the 
fractionalized Skyrmion lattice is very obvious. Natu- 
rally, the rotation can control the fractionalized Skyrmion 
lattice because it can induce the underlying vortices. 

Now, we can systematically understand the above 
questions. Its dynamics is driven by the rotation in the 
quenching process and the intrinsic spin-Hall effect de- 
rived from the effective SOC. As is well known in the 
study of rotating superfluid helium S^, [111 , the rotating 
drive pulls vortices into the rotation axis, while repulsive 
interaction tends to push them apart; this competition 
yields a vortex lattice whose vortex density depends on 
the rotation frequency. Meanwhile, the SOC causes the 
spin separation and creates the dipole structure of spin, 
which is embedded in the three- vortex structure. The 
fractionalized Skyrmion derives from the dipole struc- 
ture. In term of densities, only a core structure appears in 
the center. That is the origin of the single half Skyrmion. 
Noting that the core structure of densities agrees with 
the densities distribution of half Skyrmion in experiment 

The effect of tuning the ferromagnetic and antiferro- 
magnetic interactions. — Generally speaking, gs is much 
smaller than gn- By adjusting the two s-wave scatter- 
ing lengths ao and a2 through Feshbach resonances, the 
spin exchange interaction strength gs is tunable. We 
now perform the above experiments by changing a2. 
Figures 4(al)-(dl) show a stronger FM case of ^'^Rb 
{9s/ gn = —0.255). Contrasting the densities and the 
spin texture, we find the three- vortex structure and the 
fractionalized Skyrmion lattice are common in this FM 
BEC, though there are several half Skyrmions near the 
center. Furthermore, we also test a stronger AFM BEC 
of ^^Na {gs/gn = 0.222). The three-vortex structure 
and the fractionalized Skyrmion mainly emerge at the 
outskirt of the BECs. In fact, the strong AFM interac- 
tions restrict the fractionalized Skyrmion excitations, so 
the fractionalized Skyrmion hardly emerges in the center. 
These experiments show the universality of the fraction- 
alized Skyrmion excitation in BEC with SOC. 

Conclusion. — We have found a new type of topological 




FIG. 4: (color online). The effect of tuning FM and AFM 
interactions for BEC with k = 1 and Q = O.Scj. (al), (bl) and 
(cl) show the densities of the mp = —I, mp = I and mp = 
components of ^^Rb with ao = 101. 8as and a2 = 50.2as 
respectively (strong FM case); (dl) shows the spin texture and 
the position of vortices. The meanings of the spots and the 
colored arrows are the same as those in Fig. 2(b). (a2)-(d2) 



indicate the corresponding results of 
and a2 = llOas (strong AFM case). 



Na with ao = 50a b 



excitation-fractionalized Skyrmion-in the rotating spin- 
1 BEC with SOC. The fractionalized Skyrmion originates 
from a dipole resulted from a local spin separation. The 
fractionalized Skyrmion excitation can occur as long as 
the three-vortex structure appears, even in the strong 
FM and AFM BEC with SOC. Our study gives an ex- 
perimental protocol to observe these novel phenomena 
in future experiments. Not only do our findings exist in 
spin-1 BEC, but also the related textures should appear 
in high-spin BEC, superfluid and superconduction. This 
work is of particular significance for exploring the novel 
topological excitation such as fractionalized Skyrmions 
in quantum gas and condensed matter physics. 
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